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We consider chromatic dispersion of capillary tubes and photonic crystal fibers infiltrated with 
liquid crystals. A perturbative scheme for inclusion of material dispersion of both liquid crystal 
and the surrounding waveguide material is derived. The method is used to calculate the chromatic 
dispersion at different temperatures. 
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I. INTRODUCTION 



Together with the development of photonic crystal 
fibers (PCFs), a large amount of research has been de- 
voted to investigate the possibilities of infiltrating the 
air holes of a PCF with different liquidsQl, and thereby 
changing the optical properties of the fiber. Depending 
on the refractive index of the liquid, the guiding effect of 
the fiber can possibly be changed from guiding based 
on modified total internal reflection (mTIR), to guid- 
ing based on the photonic band gap effect, where the 
core has a lower refractive index than the effective in- 
dex of the cladding. Also selective filling of PCFs, where 
only some of the holes are infiltrated, has experienced a 
considerable interest 0, 1^ 1^, because this can be used 
to tailor the optical characteristics of the PCF. Among 
the various liquids that can be infiltrated in a PCF, liq- 
uid crystals (LC) distinguishes themselves, because of 
their anisotropic nature, which allows the possibility of 
controlling the optical parameters of the waveguide by 
changing the orientation of the moleculesQ. This orien- 
tation can be controlled in different ways, for example 
by applying an electric field externally. The optical char- 
acteristics of the fiber can also be changed by varying 
the temperature, since the ordinary and extraordinary 
refractive indices of LCs are highly dependent on tem- 
perature. Recently these tunable properties of LCs have 
been used in various experimental research projects0,0|- 
The LC can be infiltrated in the PCF using various tech- 
niques, one possibility is to use a pressure chamber, but 
this technique has shown to introduce orientational irreg- 
ularities in the alignment of LC molecules'^. Another 
possibility is to infiltrate the holes of the PCF using cap- 
illary forces, this technique has shown to give a regular 
alignment of the LC molecules. A disadvantage using 
capillary forces for the infiltration is that the length of 
the infiltrated region will only be of the order of a few cen- 
timeters, while longer infiltration lengths can be achieved 
using pressure infiltration. 



In the present work we address the problem of calculat- 
ing chromatic dispersion curves for different waveguide 
designs, where the material dispersion of both the LC 
and waveguide material is taken into account. The LC 
infiltrated PCF structures we consider have previously 
been studied theoretically jjj, without inclusion of mate- 
rial dispersion in the PCF material and LC, and only 
considering the special case where the extraordinary in- 
dex of the LC and the index of the PCF material were 
identical. The material dispersion of LC is important to 
take into account, since LCs are highly dispersive, espe- 
cially in the visible spectrum, where the dispersion can 
be much stronger than in for example silica. It has previ- 
ously been shown that approximating the total dispersion 
simply by adding the waveguide and material dispersion 
gives the correct qualitative behavior of the dispersion 
curve0, but is not sufficient if quantitative data for the 
dispersion is needed, for example to determine the posi- 
tion of zero dispersion wavelengths (ZDWs). 
To calculate precise dispersion profiles, we must there- 
fore include the material dispersion in the field equations. 
This destroys the well known scalability of Maxwells 
equations, hence if all the physical dimensions of the fiber 
are multiplied by a constant factor, we are not able to cal- 
culate the new dispersion curve without having to solve 
the field equations again. In addition, if the computa- 
tional method used to find the eigenmodes numerically 
takes the propagation constant as an input variable, and 
returns the corresponding frequency, we must ensure that 
this is done in a self-consistent manner, i.e. the values of 
the dielectric constants in the numerical calculation must 
correspond to the values of the dielectric constants at the 
frequency returned by the computational method. 
In this work we find the self-consistent frequencies based 
on a generalization of a perturbative method developed 
for isotropic waveguides |lOj. We consider dispersion pro- 
files of both simple waveguides consisting of capillary sil- 
ica tubes infiltrated with LCs, and more advanced selec- 
tively filled PCF structures. Finally we investigate how 
a change in temperature affects the dispersion character- 
istics of the fiber. 



'Electronic address: pdr@ corn^ dtu.dkl 



2 



II. THEORY 

A. Alignments of LC-molecules 

In all the fiber designs considered in this work the LC 
is contained inside a hollow circular cylinder. It is well 
known that an intense optical field will interact with the 
LC and change the orientation of the molecules [Tl|. here 
we assume that the intensity of the optical field is so 
weak that we can neglect this interaction. The LC is as- 
sumed to be in the nematic phase, where the orientation 
of the molecules is correlated, resulting in a preferred lo- 
cal orientation of the molecules. This local orientation 
is described by the director axis n(r) which is a unit 
vector pointing in the same direction as the axis of the 
LC-molecules. In the general case n(r) is found by min- 
imization of elastic energy. For the cylindrical geometry 
considered here the director axis has the following form 
n = (sin(6'), 0, cos(0)) in cylindrical coordinates (r, z). 
is the angle between the director axis n and the z-axis. 
The alignment of the LC molecules in this situation has 
previously been studied theoretically 0. 0{r) can be 
found by solving a 2nd. order nonlinear ordinary differ- 
ential equation, where the only parameters are the elastic 
constants of the LC. We assume that the molecules along 
r = are aligned parallel to the z-axis (0(0) ~ 0). The 
boundary condition at the cylinder wall depends on the 
LC and the coating of the capillary. Here we consider 
the two different possibilities 9{R) = and 0{R) — 7r/2, 
where R is the radius of the cylinder. The dielectric ten- 
sor of a nematic LC is described in terms of the perpen- 
dicular and parallel part of the optical permittivity 
and e||. If the orientation of the molecules is described 
by the angle 9{r), then the dielectric tensor has the fol- 
lowing form in cylindrical coordinates'!^ 

£00 , (1) 
e^r e^z ) 

where e^r = ej. + Ae sin^ 0, e^z = f-zr = Ae sin 9 cos 0, 
= e± and e^z = e± + Ae cos^ 9. Ae is the optical 
anisotropy defined as Ae = e|| — e^. If the 3 elastic con- 
stants of the LC describing twist, splay and bend defor- 
mations are assumed to be equal, and we further assume 
that the molecules are aligned parallel with the z-axis 
in the center of the cylinder, the orientation is given by 
9{r) =^ 2 tan^^(ar/i?), where a is a constant depending 
on the boundary condition at the wall. For the bound- 
ary conditions considered here we have the two simple 
analytical solutions 9{r) — and 9{r) = 2 tan^^(77i?), 
depending on whether the molecules are anchored paral- 
lel or perpendicular to the boundary of the cylinder. The 
two orientations are shown schematically in Fig. ^ In 
the following we will consider these two orientations of 
the LC molecules, and refer to them as planar {9{r) = 0) 
and axial {9{r) = 2 tan~^ (r/i?)) alignment. The planar 
alignment is easily achieved experimentally, while axial 
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FIG. 1: Parallel (a) and axial (b) alignment of LC 
molecules in a cylindrical geometry. Figures reproduced with 
permission!^. 

alignment requires that the capillary is coated with a 
surfactant before the LC is infiltrated|lS|. 

B. Calculation of chromatic dispersion curves 

In this section we derive a perturbative method for 
calculation of chromatic dispersion of a waveguide in- 
filtrated with LC. The method is general and can be 
applied to arbitrary waveguide designs. Our method 
is a generalization of an earlier presented method for 
isotropic waveguides |0, but this method allows the pos- 
sibility that the waveguide consists of anisotropic mate- 
rials. We consider a waveguide which is uniform along 
the z-direction, and therefore assume that the magnetic 
field can be described in the form of a monochromatic 
wave travelling along the z-direction, i.e. H(x,2;,z,i) = 
exp[i(/?z — Lot)]h.{x,y\ f}). From Maxwells equations the 
following equation for the vector field h(x, y; (3) is derived 

*h = — h, (2) 

* = V/3 X rV/jx (3) 

where the operator Vp is given by V/j = 
{d/dx,d/dy,i(3). e is the dielectric tensor, which 
in the LC region is given by the expression in Eq. JQ. 
In the silica region the dielectric tensor is simply a 
diagonal matrix, with the dielectric constant of silica in 
the diagonal. Therefore e = e{r,es{uJ),e^_{uo),e\\[lJJ)), i.e. 
the dielectric function depends on position, the dielectric 
constant of the material surrounding the LC (es), and 
the dielectric constants ej^ and e|| of the LC. Since 
material dispersion is taken into account, all 3 dielectric 
constants are assumed to be frequency dependent. The 
dispersion coefficient is defined by 

27rcWg ' ^ ^ 

where Vg is the group velocity, defined as z;g = To 
find an exact expression for the group velocity we start 
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out by rewriting Eq. |5J as 
(h,h) 



(5) 



Here and in the following wc use the following notation 
for the inner product (A,B) = / A* • Bdrj_, i.e. the 
integration is over the whole transverse plane. Now the 
group velocity is found by differentiating both sides of 
Eq. ||SJ| with respect to the propagation constant (3. The 
left hand side of Eq. (O is differentiated with respect 
to /? using the Hellman-Feynman theorem. To differen- 
tiate the operator $ with respect to /3, we note that the 
operator depends on /S explicitly through V/3, and im- 
plicitly through the dielectric constants ej — ej{uj{(3)) 
{j — S, -L,\\). In the following d/df3 denotes a differen- 
tiation for fixed dielectric constants. Using this notation 
we have the following expression for $ differentiated with 
respect to /? 



'dp 



(6) 



Where we have the following expressions for 9$ /9/3 and 



9* 

9* 

9^ 



= Vfl X e 




X e 



^V,xe-^-e- V,x 



(7) 



(8) 



Using the Hermiticity of the operator Vp, and the 
Maxwell equation V/3 x h = —iluj / c)te, now gives us 
the general expression for the group velocity in the case 
where the material dispersion e.j = ej{Lu) is known 




des 
dui 



(h,h) 

\ ' \de_L duj ^ Oeii duj ) ^ I LC 



(9) 



(10) 



here {•, •)s and {■,-)lc denote that the integration is only 
over the silica or the LC respectively. The electric field e 
is defined similarly to h, i.e. it is the part of the electric 
field where the z and t dependence has been factored 
out. In Eq. lO denotes the group velocity when the 
material dispersion is zero. An exact expression for 
is found by differentiation of Eq. lO with respect to /3, 
and again using the Hermiticity of and the Maxwell 
equation x h = —i{uj / c)ee i.e. 



2uj (h, h) 



Re{[e* X h]. 



(11) 



where (/) = J fdr±^ . Our perturbative method for calcu- 
lating D consists of several steps, first we make a guess 



for self-consistent values of the dielectric constants e^^o 
and solve Eq. From this solution we find the non- 
selfconsistent frequency wq, and the group velocity due 
to waveguide dispersion by using the definition in Eq. 
(|ll|l . Now generalizing the procedure for isotropic waveg- 
uides, we see that a first order approximation to the self- 
consistent frequency is 



Wo 



where Ej (j = 5, _L, ||) is given by 



Es 



1 (e,e)c 



2 (h,h) ' 



e)LC 



2 (h,h) • 



(12) 

(13) 
(14) 

(15) 



In Eq. (|12|l duj/dej = —^Ej has been found by differ- 
entiating Eq. with respect to ej. Acj in Eq. H12|) is 
found by noting that 



ei(^^o) + ^Uo(^sc - t^o) (16) 



= ei,o + Ae,;. 



,Y,EjAe,u;o (17) 
j 

(18) 



Eqs. H16|l define a system of 3 coupled linear algebraic 
equations. Once Eqs. Hlt)|l have been solved for Ae^, our 
approximation to the selfconsistent frequency is readily 
found using Eq. (|12|l . Since our goal is to use Eq. 10 for 
finding the group velocity, we must also find an approxi- 



mation to Vg. This is done by using that 
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A first order approximation to the self-consistent group 
velocity is then found using Eq. (|5Jl 



■ Wo 



dEj 

ap 



1 



' UJs 



1 ^3 dui 



(20) 



In Eq. (|2()|l and (|16|l we find the derivatives of the dielec- 
tric constants by differentiating the Sellmeier or Cauchy 
polynomial presented in the following with respect to fre- 
quency. Ej is found from the fields returned by the com- 
putational method, when using the dielectric constants 
Ei^o- Strictly speaking the values of Ej used in Eq. 
should be the values calculated at the selfconsistent fre- 
quency. This would require solving Eq. (|2Jl more than 
one time, and therefore significantly increase the calcu- 
lation time for each propagation constant, but we have 



i A^(25°C) 4'(25°C) A^(50°C) 4 (50°C) 



bi 



1 1.4994 


1.6933 


1.5062 


1.6395 


0.6965325 


4.368309 


10"^ 


2 0.0070 


0.0078 


0.0063 


0.0095 


0.4083099 


1.394999 


10-2 


3 0.0004 


0.0028 


0.0006 


0.0020 


0.8968766 


97.93399 





TABLE I: Parameters for Cauchy and Sellmeier polynomials 
given in Eq. I21I22II . The non-dimensionless parameters are 



all given in units of fj,m? or fj.m'^. 



found that the variations in Ej with can safely be 
neglected. The derivatives of Ej with respect to f3 are 
found using a standard three point approximation. After 
having found the selfconsistent frequencies and the cor- 
responding group velocities for a number of propagation 
constants (3, we find the dispersion D by using the defi- 
nition given in Eq. I^J, in this calculation the derivative 
of the group velocity Vg with respect to frequency uj is 
also approximated using a three point formula. 



III. RESULTS 

Both the planar and axial alignment discussed in the 
previous section are considered. For the material disper- 
sion of silica we use the Sellmeier curve 



A2 



(21) 



where ai and hi are constants. Here we use the values in 
Table m as reported by Okamoto 13]. 

For the LC a Cauchy polynomial is used for both e_L 
and e|| 
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FIG. 2: Dielectric constants of silica and E7. The curves 
for silica and E7 are based on the Sellmeier expression given 
in Eq. and the Cauchy polynomial given in Eq. H22|l 

respectively. The parameters given in Table are used. 



taken to be the values corresponding to a vacuum wave- 
length of 1/im. We solve Eq. Q using a freely available 
software package [Tsj where the electric field is expanded 
in plane waves. In this software package periodic bound- 
ary conditions are assumed on all boundaries, therefore 
all calculations for both the single capillary and the PCF 
structure are done using a supercell which is consider- 
ably larger than the LC infiltrated cylinder in order to 
minimize interactions between the images. In this work 
the distance between repeated images was 14 relative to 
the radius of the LC infiltrated cylinder for the capillary 
tubes, and 14 relative to the pitch for the PCF struc- 
tures. Each elementary cell of the supercell consisted of 
a uniform 32 x 32 grid. The relative error using these 
parameters was estimated to be below 5%, by repeating 
a set of the computations on a finer 64 x 64 grid. 



A 



-L- 



A2 



A4 



(22) 



where A^'" are constants. Here we use the values for the 
li quid crystal E7 given in Table ^ as reported by Li et 
al]^ at 25°C and 50°C. The values obtained for the di- 
electric constants using these parameters in the Cauchy 
polynomials have been shown to be consistent with mea- 
sured values throughout the visible spectrum and far into 
the infrared spectrum. In Fig. |21 we have plotted the 
three dielectric constants 65, ej^ and e|| as a function of 
vacuum wavelength. We see that the dielectric constant 
of silica [es) is below the two dielectric constants of E7 
(e^ and e||) throughout the visible spectrum and into the 
near infrared spectrum, hence a waveguide based on TIR 
can be realized in this spectrum. In the following we ex- 
amine the chromatic dispersion for different waveguides 
based on TIR for the capillary tubes, and modified TIR 
for the PCFs. Whenever the perturbative method is used 
to obtain self-consistent frequencies, the guesses for the 
self-consistent values of the dielectric constants e^-^o are 



A. Capillary tube infiltrated with LC 

First we consider a simple waveguide consisting of a cir- 
cular hole containing LC surrounded by a silica cladding. 
Such a waveguide can be realized physically by infil- 
trating a capillary tube with LC. If we assume that 
the LC molecules align in the planar orientation dis- 
cussed in the previous section, the dielectric tensor in 
Eq. only has nonzero elements in the diagonal, i.e. 
e — (im(7(e_L, e^, e||). If we further assume the cladding 
has infinite width, an analytical solution to Eq. Q can 
be derived We can therefore use this solution to 
investigate the accuracy of our numerical perturbative 
method. The fundamental mode of the waveguide con- 
sidered here is always the HEn mode. For a certain mode 
to be guided in this structure, the propagation constant 



[3 must satisfy Cg^^fc < /3 < e^'"'^, where k is the vacuum 
wavenumber k = lo/c. The fiber has a single guided mode 
when the F-parameter {V = kr{e±_ — €3102)^^'^) is less 
than 2.405, where r is the inner radius of the tube. In the 
following we consider fibers with radii of 1.5/xm, 1.0/im 



1/2, 
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FIG. 3: Chromatic dispersion of the fundamental HEn mode 
for planar alignment of LC molecules and different radii of 
the capillary tube. The lines show the chromatic disper- 
sion found using the exact analytical result for the disper- 
sion relation, and markers show the result found numer- 
ically together with our perturbative method. The inset 
shows the absolute difference between the two results, i.e. 

l^-^-^l — \D exact Dpej^f_urhative\- 



and 0.75/im, these fibers are single mode for wavelengths 
larger than 1.63/im, 1.06/iin and 0.81/im respectively. In 
Fig. 13 we have compared the chromatic dispersion found 
analytically with the chromatic dispersion found using 
the numerical method described above, together with the 
perturbative method described in the theory section. We 
see that the dispersion curves found numerically together 
with the perturbative method are quantitatively consis- 
tent with the exact dispersion curves. In the inset in Fig. 
Othe difference between the analytical result and the per- 
turbative result is also plotted, we see that the smallest 
deviations between the two results occur in the infrared 
region, this is also expected since the material dispersion 
is lowest in this region (see Fig. [^J . The relative error of 
the perturbative method is below 5% in the wavelength 
interval from 0.6/im to 1.7/im. The dispersion is very 
high in the visible spectrum, which is mainly due to the 
high material dispersion in this region. Also notice that 
the dispersion is normal {D < 0) for all the wavelengths 
and radii considered for the planar alignment. 
For the axial alignment of the LC molecules there only 
exists an analytical solution to Eq. Q for the TE 
modesp^. But since the fundamental mode, i.e. the 
mode with the lowest frequency, is not a TE-mode we 
must solve Eq. (0) numerically to find the chromatic dis- 
persion for the fundamental mode. It turns out that at 
short wavelengths the fundamental mode is the TMqi 
mode, while for longer wavelengths the fundamental 
mode is the hybrid HEn mode. For the tube radii and 
wavelengths considered here, the capillary tubes with the 
axial orientation are always multimoded. Here we con- 
sider the two modes with lowest frequency; the HEn and 
TMqi mode. The dispersion of these modes as a function 
of vacuum wavelength is shown in Fig. ^ Again we see 
that the dispersion is mostly normal for the wavelengths 
and tube radii considered here. But for r = 1.5/im the 
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FIG. 4: Chromatic dispersion for axial alignment of LC 
molecules and different radii of the capillary tube. The 
self-consistent frequencies are found using the perturbative 
method. Examples of |Hp for the TMqi and HEn mode are 
shown on the insets. 




FIG. 5: Representative contour plots of real part of the x- 
component of H for the TMqi mode (left) and the HEn mode 
(right) for the axial alignment in a capillary tube. The white 
circle shows the boundary of the LC infiltrated capillary. 



dispersion becomes anomalous for vacuum wavelengths 
higher than approximately 1.4/im for the TMqi mode. 
For r = 1.0/im the fundamental mode switches from the 
HEii mode to the TEoi mode at a vacuum wavelength 
around A = 1/im. For r — 0.75^m and r = 1.5/im the 
switch between the two modes happens below A = 0.6/im 
and above A = 1.7/im respectively. In an experimental 
setup light is coupled into the LC infiltrated region using 
the HEii mode of a single mode step index fiber which is 
an even mode. Therefore it will most likely be easiest to 
excite the HEn mode of the LC infiltrated region, since 
this mode is also even, in contrast to the TMqi which is 
odd. This is demonstrated in Fig. [S] where the real part 
of the x-component of H is plotted for the TMqi and 
HEii. A similar behavior is found for the other compo- 
nents of the H-field, hence the TMqi mode is odd (even 
though the intensity plot of |Hp in Fig. 21 is even), and 
the HEii mode is even. 

In the following the effect of increasing the temperature 
to 50°C will be studied. We do not consider tempera- 
tures above 50°C, since the clearing temperature, i.e. the 
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FIG. 6; Chromatic dispersion at 25°C (solid lines) and at 
50°C (dashed lines) for HEn modes, (a)-(b) Capillary tube 
{R = 1.0/im) with LC parallel and axially aligned respectively. 
(c)-(d) Selectively filled PCF structure as shown in inset in 
Fig. |7|with LC parallel and axially aligned respectively. In 
subplot (c) the bullets (•) and circles (o) are for R = 0.75/im 
and R = 1.0/im respectively. The dispersion profile in (d) is 
for R — 1.0 fim. 



temperature where e±_ = e||, is around 57°C for E7|l4|. 
Above the clearing temperature the LC is no longer in 
the anisotropic nematic phase. The parameters for the 
ordinary and extraordinary indices of refraction at 50°C 
can also be found with Cauchy polynomials. The coef- 
ficients at 50°C are given in Table HJ and the dielectric 
constants at 50°C are plotted in Fig. |21 We see that the 
increase in temperature also increases the ordinary dielec- 
tric constant, while the extraordinary dielectric constant 
is lowered. In subplot (a) and (b) in Fig. Elthe effect of 
raising the temperature to 50°C is shown for the capil- 
lary tube with the planar and axial alignment of the LC 
molecules. We see that the dispersion increases for both 
alignments. For the planar alignment the explanation 
for this is straight forward. The HEn mode carries most 
of its energy in the transverse components of the field, 
and since the transverse components experiences the or- 
dinary dielectric constant, which increases with temper- 
ature, the temperature increase effectively increases the 
index difference between the core and cladding. The in- 
creased index difference gives rise to the higher disper- 
sion. For the axial alignment the explanation for the in- 
creased dispersion is more complicated than for the pla- 
nar alignment. Here a field which is mostly transverse 
will experience e± near the center of the cylinder, and e|| 
near the wall of the cylinder. Since e± increases and e|| 
decreases with temperature, as shown in Fig. |21 it is dif- 
ficult to say a priori whether the dispersion is increased 
or decreased. 



B. PCF infiltrated with LC 

In this section we consider a PCF-design similar to the 
structure recently investigated by Zografopoulos et a/Q. 
Where the possibility of changing the fiber characteristics 



FIG. 7; Chromatic dispersion for the PCF structure shown 
in the inset with difi'erent radii of the central hole. The ratios 
-R/A = 0.6 and r/A = 0.3 are fixed. The LC molecules are 
planarly aligned. A zoom around the ZDWs is shown in Fig. 

m 



-R=1.5nm (TM) 
' R=1.5nm (HE) 
-R=1.0jim (TM) 
I R=1.0nm (HE) 
-R=0.75nm (TM) 
-R=0.75nm (HE) 




1.6 0.8 1 1.2 1.4 1. 
X [|im] 



FIG. 8; Chromatic dispersion for the PCF structure shown 
in inset in Fig. |7|with difi'erent radii of the central hole. The 
LC molecules are aligned axially. Examples of |H|^ for the 
TM and HE mode are shown on the insets. 



by applying an external electric field was considered. The 
structure has a cladding consisting of airholes placed in a 
triangular structure and a core hole infiltrated with LC. 
A cross section of the considered structure is shown in the 
inset in Fig. A physical realization of such a struc- 
ture will require selective filling, which has recently been 
demonstrated 0, Q, where selective filling was achieved 
by collapsing the small holes using a fusion splicer, and 
then the holes with the larger radius were infiltrated. 
Here we consider a structure where the infiltrated cen- 
ter hole has a radius twice as large as the radius of the 
cladding holes. The pitch is 5/3 times the radius of the 
central hole. Compared to the capillary tube studied 
in the previous section this structure has a higher index 
contrast between the core and cladding, because the pres- 
ence of the airholes significantly lowers the effective index 
of the cladding. Again we consider both the planar and 
axial orientation of the LC in the center hole. In Fig. 
13 the dispersion curves for different radii of the center 
hole are shown for the fundamental mode HE-mode with 
the planar alignment of the LC. The fiber is multimoded 
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for the wavelengths considered here. We see that all the 
fibers now have regions of both normal and anomalous 
dispersion. The fiber with a center hole radius of 0.75/zm 
has two ZDWs at A = 1.05/im and A = 1.55/im. The 
fibers with center hole radii of 1.0/im and 1.5/im each 
have one ZDW at A = 1.15/im and A = l.Sfim respec- 
tively. The dispersion curves for the axial alignment of 
the LC are shown in Fig. |H1 Like the capillary tube with 
the axial alignment, the type of the fundamental mode 
is also dependent on the wavelength for the PCF with 
the LC axially aligned. We see that the fibers with cen- 
ter hole radius 1.0/im and 1.5/im, now have large regions 
where the dispersion is anomalous for the mode that re- 
sembles the TMoi mode of the single capillary tube. For 
the mode that resembles the HEn-mode of the single 
capillary tube the dispersion is purely normal for all the 
waveguide designs considered here. 

The effect of increasing the temperature is also investi- 
gated for this waveguide design. In plots (c) and (d) in 
Fig. El we see again that the dispersion increases with 
temperature. In subplot (c) the dispersion profiles for the 
PCFs with R = 0.75/im and R = 1.0/im are shown. The 
lower ZDW of the fiber with R = 0.75/im can be tuned 
between 1.051/im at 25°C and 1.058/im at 50°C, while 
the higher ZDW can be tuned between 1.565/im and 
1.571/im. The fiber with R — 1.0/im has one ZDW in the 
optical spectrum which can be tuned between 1.128/im 
and 1.139/im. For the PCFs with the axial alignment of 
the LC molecules we do not have anomalous dispersion 
for the designs considered here. But the plots for for the 



axial alignment in Fig. indicates that the dispersion 
can be tuned in a broader interval for this alignment. 



IV. CONCLUSION 

An accurate method for calculating chromatic disper- 
sion of anisotropic waveguides is demonstrated. The 
method is based on a generalization of a previously pre- 
sented method for isotropic waveguides. We have ap- 
plied the method to a simple step index fiber with an 
anisotropic LC core since this problem has an analyti- 
cal solution. Our results show that the method can be 
applied to calculate chromatic dispersion curves that are 
consistent with the exact result throughout the visible 
spectrum and into the near infrared spectrum. 
With the method we have studied chromatic dispersion 
of capillary tubes and PCFs infiltrated with LC. The con- 
sidered PCFs are all multimoded in the wavelength in- 
tervals considered, while it is shown that single mode op- 
eration is possible for the capillary tube infiltrated with 
LC molecules aligned in parallel. The tunability of the 
different LC infiltrated waveguides is investigated by cal- 
culating the chromatic dispersion at 25°C and at 50°C. 
For the two different alignments of the LC considered 
here, the tunability is highest for the axial orientation, 
while the tunability for the planar orientation is weaker. 
A waveguide design where the ZDWs can be tuned over 
approximately lO/xm is demonstrated. 
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